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The superfluid Stirling refrigerator (SSR) is a Stirling machine that uses the 
ideal gas behavior of the 3 He Component of a 3 He- 4 He mixture to provide 
cooling at temperatures below 2 K. The design of SSR’s to date has used the 
Schmidt model to predict the performance of a given machine. Unfortunately , 
the Schmidt model does not account for the phonon and roton excitations 
present in 4 He component of the mixture. These excitations significantly 
change the performance of the SSR. An analytical model of the SSR that 
accounts for the phonon and roton excitations is presented. High temperature 
SSR experimental data is also presented and found to compare well with the 
model. 


INTRODUCTION 

Superfluid Stirling refrigerators (SSR) provide cooling below 2 K, by 
exploiting the unusual properties of superfluid 3 He- 4 He mixtures. 1 To first 
approximation, the 3 He component of the mixture behaves as an ideal 
Boltzmann gas in an inert background of superfluid 4 He. By using super¬ 
leak bypassed pistons, the 3 He component of the mixture is driven through 
the Stirling cycle to provide cooling. SSR’s have demonstrated cooling to 
0.17 K 2 and have demonstrated cooling powers of 977 at 0.5 K. 3 

The zeroth order design of all Stirling machines is done with the 
highly idealized Schmidt model, 4 which assumes an ideal Boltzmann gas 
behavior for the working fluid. This model is extremely useful as a bench¬ 
mark for all Stirling machine designs and provides a framework for under¬ 
standing the overall trends in the design space. Unfortunately, the ideal 
Boltzmann gas model, and consequently the Schmidt model, is appropriate 
only for a limited temperature and concentration range of 3 He in 4 He 
which, in turn, makes the Schmidt model of limited use for superfluid 
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Stirling refrigerator design. At high temperatures, the thermal excitations in 
the 4 He, the phonons and rotons, significantly alter the thermodynamic 
and kinematic behavior of the working fluid from that of an ideal 
Boltzmann gas; whereas at low temperatures, the 4 He excitations are 
absent but the 3 He component of the mixture behaves as a Fermi, not a 
Boltzmann, gas. 

Brisson and Swift 5 have shown that the presence of these phonon- 
roton excitations in the 4 He component of the mixture can significantly 
augment the cooling power of the SSR. However their analysis and data 
were limited to SSR operation at uniform temperature. Here, we will 
extend the results of Brisson and Swift to the non-uniform temperature 
operation of an SSR where the phonon-roton effects are incorporated into 
an Schmidt-like analytical model. We will not, however, address in detail 
the low-temperature Fermi behavior of the 3 He here. 

We believe that the phonon-roton (PR) model presented here is the 
appropriate zeroth order model for superfluid Stirling refrigerators. 
Because the model is of closed form, it is easily programmed onto a 
spreadsheet for design evaluation. In addition, it predicts several effects not 
predicted by the Schmidt model, such as a non-zero ultimate temperature 
and an enhanced cooling power above 1 Kelvin. 

The paper is divided into several sections. First, a simple discussion of 
the basic Stirling cycle is described and how it relates to actual superfluid 
Stirling refrigerator designs. This is followed by a review of the Schmidt 
model for Stirling machines that establishes notation and provides a 
framework for the subsequent description of the phonon-roton model. The 
experimental SSR and the experimental procedure are then described and 
the experimental results are compared to the Schmidt and PR models. The 
penultimate section of this paper describes other limits to SSR performance 
such as heat flush, Fermi, and phase separation effects. The paper ends 
with some closing comments. 


Cycle Basics 

The superfluid Stirling cycle can be understood using Fig. 1, which 
depicts the four steps in an articulated Stirling cycle. The liquid 3 He- 4 He 
mixture is held between a cold piston and hot piston, each bypassed by 
superleaks. The volume behind of each of these pistons is filled with liquid 
superfluid 4 He. The piston cylinders are connected by a regenerator which 
can be thought of as an array of high heat capacity tubes with poor ther¬ 
mal conductivity in the longitudinal direction and good thermal conduc¬ 
tivity in the transverse direction. (In the ideal regenerator the heat capacity 
of the tubes is infinite, the thermal contact between the fluid and the tubes 
is perfect and there is no thermal conductance in the longitudinal direction.) 
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Fig. 1. A schematic diagram of the four stages of an articulated superfluid 
Stirling cycle. The compressor (expander) piston is assumed to be isothermal and 
in good thermal contact with a heat reservoir at high (low) temperature T, (T,,). 


We assume for the moment that the fluid volume contained in the regener¬ 
ator is negligible when compared to the fluid in the piston cylinders. 

The cycle begins with the 3 He- 4 He mixture entirely in the compressor 
piston cylinder; at this point, the mixture is at the same temperature as the 
high temperature reservoir, T,. The normal 3 He component is isothermally 
compressed by the hot piston (Fig. la) while the superfluid 4 He component 
of the mixture flows freely through the superleak to the backside of the 
compressor piston. Because the 3 He component of the mixture behaves as 
an ideal gas in an inert “atmosphere” of 4 He, heat is rejected to the high 
temperature reservoir (— Q , in Fig. la). Following this, both pistons are 
slowly moved in tandem to isochorically displace the 3 He through the 
regenerator to the cold expander piston cylinder (Fig. lb). Ideally, there is 
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a linear temperature gradient in the regenerator and, because of the large 
heat capacity of the regenerator matrix, the 3 He is reversibly cooled from 
the compressor piston temperature to the temperature of the cold thermal 
reservoir, T e . The 3 He is now isothermally expanded, absorbing heat, Q e , 
from the low temperature reservoir (Fig. lc). Finally, the expanded 3 He is 
displaced back through the regenerator (where it reversibly reclaims the 
heat deposited there earlier) to the compressor piston cylinder (Fig. Id). 
The refrigerator is back in the original state and the net effect has been a 
heat transfer from the cold reservoir to the hot reservoir. 

A schematic diagram of a more practical SSR, first demonstrated by 
Brisson and Swift, 6 is shown in Fig. 2. The refrigerator is actually two 
Stirling refrigerators operating 180° out of phase with each other. The com¬ 
pressor pistons are driven sinusoidally, 90° out of phase with the expander 
pistons. This back-to-back Stirling configuration allows the use of a 
recuperator (a counterflow heat exchanger) instead of a regenerator to 
provide the thermal isolation between the compressor and expander piston 
cylinders. The primary reason for a recuperative design over a regenerative 
design is the dearth of sub-Kelvin high-heat-capacity materials necessary 
for an efficient regenerator. This design also allows the removal of the pure 
4 He space shown in Fig. 1 since each Stirling refrigerator half acts as a 
reservoir of 4 He for the other refrigerator. 

The modeling of a recuperative Stirling machine (Fig. 2) can be 
reduced to a model of a regenerative Stirling machine (Fig. 1) by modeling 



piston-cylinders piston-cylinders 

(high temperature) (low temperature) 

Fig. 2. A schematic diagram of a practical superfiuid Stirling cycle. The mechanical 
linkage between the compressor pistons insures 180° phase shift between the com¬ 
pressor pistons. The expander pistons are similarly configured. The compressor and 
expander pistons are driven sinusoidally and 90° out of phase with each other. To first 
order, the flows in the recuperator are equal and opposite as shown. 
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each half of the recuperative Stirling machine as a regenerative machine. 
The flow on one side of the recuperator acts as an infinite heat capacity for 
the flow on the other side of the recuperator. So, from the perspective of 
the flow on one side of the recuperator, the recuperator looks like a 
regenerator. Since both the Schmidt model and the phonon-roton model, 
model only one side of the recuperative Stirling machine we will refer to the 
recuperator as the regenerator in these sections. 


Schmidt Model of the SSR 

We review here the Schmidt model for the superfluid Stirling refrigerator. 
This model is the zeroth order model for the design of any Stirling machine 
and was used extensively in the original designs of SSR’s. The Schmidt model 
is an ideal model; there are no viscous effects and all heat transfers occur 
reversibly. Consequently, the efficiency of this cycle is that of Carnot. We out¬ 
line the model here for completeness, to allow its comparison to the phonon- 
roton model, and to provide a simple framework for the understanding of the 
more complicated PR model. This treatment borrows liberally from Ureili’s 
discussion of the Schmidt analysis of Stirling machines. 7 

Here, we analyze one side of the back to back superfluid Stirling 
refrigerators. This is entirely equivalent to modeling the SSR depicted in 
Fig. 1. The temperature of the liquid throughout each piston-cylinder is 
assumed to be at the temperature of the local heat reservoir. We also 
assume that the 3 He behaves as an ideal Boltzmann gas in an inert back¬ 
ground of superfluid 4 He. The volumes of the piston cylinders are assumed 
to be driven sinusoidally, namely as 


V = V, + V 

r c r etc ' ' swe 


1 -l- cos 9 


V - V, + V 

r e r cle 1 swe 


1 + cos (6 + a' 


( 1 ) 

( 2 ) 


where V c , V clc , V swc are the compressor piston’s total volume, clearance 
volume and swept volumes, respectively. Similarly, V e , V c/e , V swe are the 
expander piston’s total volume, clearance volume and swept volume, 
respectively. The crank angle, 6, is assumed to vary linearly with time, 
6 = a>t, where co is the angular frequency and t is time. There is a phase 
shift between the pistons, a, of 90°. 

Because of the infinite heat capacity of the ideal regenerator matrix 
and, by assumption, the regenerator is in steady state, the temperature dis¬ 
tribution in the ideal regenerator is independent of time. Since the mass 
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flow through the regenerator is balanced (the net mass flow through the 
regenerator is zero) and the heat capacity of the ideal gas is a constant, the 
temperature in the regenerator varies linearly with position from the com¬ 
pressor temperature to the expander temperature. 

Since the pressure drops associated with the viscous flow of the ideal 
gas are assumed negligible, the pressure is spatially uniform throughout the 
Stirling machine but time varying. The total number of moles of 3 He, N 3 , 
in the Stirling refrigerator can be written as: 




RT^ 


rt e 


(3) 


where p is the time varying pressure, R is the universal gas constant, V r is 
the regenerator volume and, T c and T e are the compressor and expander 
operating temperatures. T LM is the log mean temperature of the 
regenerator (=(T C — T e )/\n(T c /T e )). From Eq. (3), the pressure associated 
with the 3 He is then 


P = 



(4) 


The time dependence of p is buried in the time dependent variation of V c 
and V e . The work per cycle done by the expander and compressor is deter¬ 
mined by integrating pdv for each piston over one cycle. The results are: 


where 


W r = nV n 


W r = n V„ 


NR 


W 1 ~ b s 
NR 


(^) 


sin(/0 


— a) 


(5) 

( 6 ) 


V, V 

g _ ClC swe 

s ~ T r 2T„ 


V r 

+ ~+~+ 


c/e K, 


2 T' 


P = tan 1 


V swe sin cl/T € 


V swe cos a/T e + V swc /T t 


1 

C ~ 2' 


T, 


2V V 

, swe swe _ . 

+ ^x““ s “ + 


V 

swe 


(7) 

( 8 ) 

(9) 



(10) 
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The heat transfers to the gases in the expander can be calculated by 
applying energy conservation to the expander. For an ideal gas, the energy 
equation is 

<^ = Q-W+N iCp T ( 11 ) 


where E is the internal energy of the 3 He in the control volume, Q and W 
are the rate of heat and work transfer, positive N 3 is the molar flow rate 
of 3 He into the expander volume, c p is the molar heat capacity at constant 
pressure (=5/2 R) of the 3 He, and T is the temperature of the 3 He entering 
and exiting the expander. As we are dealing with ideal components, the 
temperature of the gas entering and leaving the expander is always T e . 
If Eq. (11) is integrated over one cycle, the N 3 c p T term reduces to a con¬ 
stant, c p T e , times a cyclic integral of the molar flow into and out of the 
control volume. This integral is zero by conservation of mass. Similarly, the 
internal energy term integrates to zero because the starting and finishing 
states are the same. We conclude that the heat absorbed by the expander 
is equal to the work done by the expander or Q e = W e . A similar analysis 
on the compressor results in the conclusion that the heat absorbed by the 
compressor is equal to the work done by the compressor or Q c = W c . 

The coefficient of performance of the refrigerator, COP, is defined as 
the heat transferred to the expander divided by the total work absorbed by 
the refrigerator. Using Eqs. (5) and (6), we find the coefficient of perfor¬ 
mance of a Schmidt cycle refrigerator is 


COPs - 


Qe 

w c + w e 



which is the Carnot coefficient of performance. 


( 12 ) 


Phonon-Roton Model 

Introduction 

In the discussion of the Schmidt model for the SSR we blithely assumed 
that the contribution of the 4 He component of the mixture does not con¬ 
tribute to the performance of the SSR. Here, we account for the thermo¬ 
dynamic contribution of the normal component of the 4 He. Our goal is to 
determine the work and heat transfer in closed form for each of the SSR 
piston cylinders. To do this, we follow the same overall procedure described 
in the Schmidt model section. First, an expression for the pressure of the 3 He 
in the SSR is developed. The effect of the phonon-roton gas on the 3 He 
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concentration in the SSR is considered and an expression for the effective 
pressure on the SSR pistons is determined. This is compared to the corre¬ 
sponding expressions in the Schmidt model. Following this, the expressions 
for the work done on the pistons are determined by integrating the pressure- 
volume relations. To evaluate the heat transfer of the pistons, the first law is 
introduced. But, unlike the Schmidt model, there is non-zero contribution to 
the energy equation due to a net (osmotic) enthalpy flow through the regen¬ 
erator. We then develop a simple model for the regenerator that allows the 
evaluation of the enthalpy term in the expander piston and subsequently the 
heat transfer in the expander piston-cylinder. The heat transfer in the com¬ 
pressor piston is then determined by applying the first law to the entire SSR. 
Finally, since the entropy and fountain pressure of pure 4 He are needed to 
evaluate the derived heat and work transfer expressions, we develop an 
empirical expression for these quantities by fitting to experimental data. 

3 He Pressure in the SSR 

If critical velocities are not exceeded, the superflows in the 4 He compo¬ 
nent will maintain a constant 4 He chemical potential throughout the SSR. 
Barring viscous effects, force balance requires that the pressure be uniform 
throughout the 3 He- 4 He mixture. Across the superleaks, the 4 He chemical 
potential is constant. If the mixture is in contact with a large reservoir of pure 
4 He at pressure P r and compressor temperature T c , the 3 He component will 
distribute itself to satisfy the thermodynamic relation 

11° 4 {P„ T c )=fi 4 {P r + n, T c , x c )=n 4 (P r + n, T e ,x e ) ( 13 ) 

where /,i° 4 and /i 4 are the 4 He chemical potential in the pure state and in the 
3 He- 4 He mixture, respectively, n is the osmotic pressure, T e is the tem¬ 
perature of the expander cylinder, and x c (x e ) is the molar concentration of 
the 3 He in the compressor (expander) cylinder. In the actual SSR this large 
reservoir of pure 4 He is fictional; however, the effect of the opposing SSR is 
as if a large volume of pure 4 He exists at the center of the superleak between 
the two SSRs. The first equality of Eq. (13) is the defining relation 
for the osmotic pressure. The second equality relates the molar concentra¬ 
tions found in the compressor and expander cylinders for a given T c and T e . 

We will assume here, contrary to the Schmidt model, that the regener¬ 
ator volume is negligible. This simplifying assumption will be discussed in 
more detail later. The conservation of mass (moles) of 3 FIe in the space 
enclosed by the superleaks results in 

N 3 = f — dV+ f — dV=(— f) +(— f) (14) 

J c v m 4 v m \v m j c \v m J e 
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where N 3 is the total number of moles of 3 He in the refrigerator, x is the 
3 He molar concentration, v m is the molar volume of the mixture, and the 
subscripts c and e refer to evaluating the associated quantities with com¬ 
pressor or expander values. Equations (13) and (14) determine the 3 He 
concentrations in the refrigerator. To proceed further we require a relation 
for the 4 He chemical potential of the mixture. Kuerten and co-workers 8 
suggest that the 4 He chemical potential of a mixture can be written as 

» 4 (p,T,x)=»° 4 (p,T)-n(T,x)v ° 4 (15) 

where is the molar volume of pure 4 He and p is the pressure of the 
mixture. The latter part of Eq. (13) can be expressed as the differential of 
p 4 is zero throughout the 3 He- 4 He space or using Eq. (15): 

(tp 0 4 -v°dn = v°dp-s 0 4 dT-v°dn= - s° 4 dT - v °dU = 0 (16) 

where 5 “ is the molar entropy of pure 4 He. The pressure term, vj dp , is zero 
since mechanical equilbrium requires that there are no pressure gradients 
in the 3 He- 4 He part of the device. The expression can now be integrated to 

[ P 4 .S 4 dT T n(x, T) = constant (17) 

J o 


where is the molar density of pure 4 He, The first term is the fountain 
pressure of pure 4 He which we will label as p f . The osmotic pressure is 
identified as the pressure of the quasiparticle gas p q = Tl(x, T). 9 Equation 
(17) reduces to 


(pAT)+P q (x,T)) c = (p f (T)+p q (x, T)) e . (18) 

Previous work on the SSR has always arranged to make the p f terms 
in Eq. (18) insignificant. This is true at low temperatures for reasonable 
concentrations (T< 1 K, x~0.06) because the fountain pressures are much 
smaller than the quasiparticle gas pressures. Watanabe et a /. 10 made high- 
temperature measurements that eliminate the mechanical fountain effects 
by running the SSR at uniform temperature since, in this case, the fountain 
terms cancel out identically. 

Since our intent is to develop a simple Schmidt-like model for the high 
temperature operation of the SSR, we will assume that the 3 He quasipar¬ 
ticle gas behaves as an ideal Boltzmann gas, p g = n 3 RT where n 3 is the 3 He 
molar density in the mixture. The fountain pressure can be thought of as 
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the pressure of a phonon-roton gas of quasiparticles present in the 4 He 
component of the mixture. Equation (18) becomes 

n 3e RT e -n 3c RT c =p f (T c )-p f (T e ). (19) 

(The Schmidt model assumes that the right hand side of Eq. (19) is negli¬ 
gible and proceeds from this point.) If N 3 moles of 3 He are loaded into the 
SSR then Eq. (14) becomes 


N 3 = n ic V c + n 3e V e . (20) 

Equations (19) and (20), with the ideal gas relation for p q , can be 
manipulated for the quasiparticle gas pressures in the compressor and the 
expander volumes, resulting in: 


/ > v * 
N 3 R [Pfc Pfe) T e 

P9°— y v ~ V V 

I c . Y e 

T ^ T T ^ T 

1 c 1 e x c x e 


Pqe 



v 

( Pfe-Pfe)TjT 
x c 



( 21 ) 


( 22 ) 


where p fc =P/i T c ) and p fc = p f { T e ). The first terms in these equations 
would be the result if the fountain pressures were ignored and corresponds 
to the Schmidt pressure of Eq. (4). In each of these equations, the second 
term is a modification in the pressure due to the fountain effect flushing 
3 He atoms from the compressor piston cylinder into the expander piston 
cylinder (assuming T c > T e ). Since the piston motions for the PR model are 
identical to the Schmidt model, Eqs. (1) and (2) can be substituted for the 
volumes in the numerators of Eqs. (21) and (22) resulting in 


Pqc 


P 


qe 


N 3 R 


A P 

T, 


fly , ’ swe 


v c V e 
— + — 
T„ T„ 


N 3 R+^[v clc +^ 
y y 

_£ i _£ 

T T 


AP f V swe cos(6 + <x) 



AP/V mc cos(B) 
2 T c Vc V. 
T c + T e 


(23) 


2 


(24) 
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where A P f is the fountain pressure difference between the compressor and 
expander pistons (A P f =p fc —p fe ). There is a time dependence in all the 
terms of Eqs. (23) and (24) (through 8 , V c (0) and V e {0))\ however, the 
numerator of the first terms is a constant in time. These first terms are 
pressure oscillations that have a Schmidt character to them. The com¬ 
pressor piston has a reduced effective 3 He density whereas the expander 
piston has an enhanced effective 3 He density due to the fountain pressure. 
The second terms correspond to a modulated pressure oscillation directly 
related to the fountain pressure difference between the pistons. Equations 
(1) and (2) can be substituted for the volumes in Eqs. (23) and (24) result¬ 
ing in 

(^R)eff.c AP f V swe cos(8 + a) 

P qc .v(l+/i cos <p) 2 T e s(\+b cos tj>) 

(A 'R)eff,e A P f V swc COS (8) 

qe s(\+b cos (j>) 2 T c ,v( 1 + b cos (f >) 

where 


(NR) e// , c = N 3 R-^( V cle + ^ 
(NR) e/f , e = N 3 R + ^^V clc + ^ 



</>^0 + p 
b = c/s 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


The substitutions outlined in Eqs. (25)—(33) reduce the pressures oscilla¬ 
tions to simple form where the time dependence is entirely contained in the 
arguments of the trigonometric functions through <j) and 0. We have 
repeated the definitions of /f and c for purposes of clarity. 
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Work Transfer to the Pistons 

The work done by the gas can be calculated for each piston by directly 
integrating pdv over the entire cycle. For the compressor piston, the work, 
W c , is 

IV --( 2n ( (NR> ^c _ APf Vswe cos(fl + «) \/ V swc sin 6 \ 

c Jo \s(l + b cos </>) 2T e s(\+b cost) A 2 ) 1 J 


Similarly, the expander piston work, W e , is 


f 2 * / {NR) eme AP f V swc cos(0) V V swe sin(6> + a) 

o \s(\ + b cos <j>) 2T C s(\ + b cos <f>))\ 2 


dO. (35) 


The integrals can be evaluated by making a substitution of y = tan(#/2) 
and integrating by partial fractions. The results are: 


nV 

W,z= swc 




i w ( NR) e/ y c . AP fV s . ,^ n ^ 

1 If- ])~ sm (P) ~ 2 Tb 2 sin(2 ^~ a) 


+ - 


nV„ 


. cos(/?) sin(/? — a) 


iT'jr^b 2 


(36) 


and 


7T V 

w - “ Lm 


* V n 


^ sinW _ a)+ ^^£ Sin(2/; _ a) 


APj ~ sw c —.s\n(p) cos(fj — a.) 


2 r c yr^b~ 2 


(37) 


Conservation of Energy Applied to the Expander 

We have determined the work done by the expander and the com¬ 
pressor pistons. In the case of the Schmidt model, we argued that these 
works are the same as the heat absorbed by each of the pistons. Unfor¬ 
tunately, this is not the case with the presence of the phonon-roton gas. 
Here, we evaluate the heat transfer experienced by the expander piston by 
applying the first law to the expander piston. 

Ebner and Edwards 11 have shown that the conservation of energy 
equation for steady-flow no-shaft-work systems in helium mixtures can be 
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cast in terms of the osmotic enthalpy. Their arguments are easily extended 
to the unsteady-flow with shaft-work case. The result is 

dE 

= Q-W+ £ N 3 h? + X (38) 

ports ports 

where E is the total energy in the control volume, Q is the rate of heat 
transfer into the control volume, W is the rate of work transfer out of the 
control volume, N 3 is the 3 He molar flow rate into the control volume, 
h° s is the osmotic enthalpy per mole of 3 He, 7V 4 is the molar flow rate of 
the 4 He into the control volume and p 4 is the 4 He chemical potential. The 
indicated sums are over the ports that allow the mixture to flow into and 
out of the control volume. We have assumed in writing Eq. (38) that the 
gravitational and kinetic energies are not significant. Applying Eq. (38) to 
the expander piston cylinder and integrating over a complete cycle, the 
conservation of energy reduces to 

Q,= W.-jhfdN 3 (39) 

where Q e is the heat absorbed per cycle by the expander, W e is the work 
per cycle done by the expander. The dE/dt term in Eq. (38) does not con¬ 
tribute because the internal energy is the same at the beginning and end of 
the cycle. The X N 4 p 4 term also does not contribute to the integral. The 
4 He chemical potential is held constant by the large 4 He reservoir in con¬ 
tact with the SSR through the superleaks. Since the 4 He chemical potential 
is constant throughout the cycle, the cycle integral reduces to summing of 
all the 4 He molar flow into and out of the control volume over a complete 
cycle. This is zero by conservation of mass. 

In order to find Q e , it is necessary to evaluate the cyclic integral of the 
osmotic enthalpy, which is defined as 11 

h 3 =Hi +~ = M q + (P~P q ) v 3 + — = h q + [p-p q ) V 3 T~—— Ts° 

X X X 

(40) 

where p 3 is the 3 He chemical potential, T is the temperature, s m is the 
molar entropy of the mixture, p q is the chemical potential of the 3 He 
quasiparticle gas, p is the pressure of the system, r 3 is the partial volume 
of the 3 He, x is the molar concentration of 3 He, h q is the molar quasi¬ 
particle gas enthalpy. The cyclic osmotic-enthalpy integral in Eq. (39) 
reduces to three integrals: an integral of the 3 He quasiparticle gas enthalpy 
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(§hgdN 3 ), a flow work integral ($ {p —p q ) v dN 3 ), and an entropy flow 
integral (f ((1 —x)/x) Ts^dN 3 ). 

It is easily shown that the flow work integral has value zero. The par¬ 
tial volume of the 3 He is related to the concentration dependence of the 
molar volume of the 3 He- 4 He mixture. DeWaele 12 suggests that the molar 
volume of the mixture, v m , can be written as: 

v m = v°(l +a'jc) (41) 

where is the molar volume of pure 4 He (=27.58 x 10“ 6 m 3 /mol) and 
a'= 0.286. The partial volume is calculated using the partial volume rela¬ 
tion 13 for binary mixtures: 


v 3 = v m + (\-x)-^=v%(l+a'). (42) 

We see then that the partial volume of 3 He is constant. We have assumed 
that the SSR is connected to a large volume of 4 He at constant pressure, 
P r , and temperature, T r . As suggested by Eq. (13), the difference p—p q is 
equal to the constant P r . The cycle integral in Eq. (39) of the (p—p q ) v 3 
term reduces to a constant times the cycle integral of the number of 3 He 
particles in the volume and by conservation of mass, the flow work integral 
is zero. 

The other two integrals cannot be dispatched with the same ease as 
the flow work integral because these integrals depend on the temperature 
of the fluid entering and leaving the expander. In order to evaluate these 
integrals, we will find in the next section that the temperature of the fluids 
entering and exiting the expander is at the time invariant temperature of 
the expander’s thermal reservoir ( T e ). 


The Regenerator 

In modeling the regenerator in the above cycle we have assumed that 
the regenerator has zero volume and is 100% effective. This is actually 
more restrictive than the standard Schmidt model, which allows a non-zero 
volume in the regenerator and compensates for this volume in the conser¬ 
vation of mass equation by using a log-mean temperature for the 
regenerator temperature. Unfortunately, due to the thermal imbalance 
introduced by the phonons and rotons present in the fluid mixture, the 
temperature distribution in the recuperator of the SSR is spatially non¬ 
linear. A thorough treatment of the regenerator is beyond the scope of this 
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work; however, several simplifying assumptions can be used in our 
phonon-roton model. 

Due to the presence of the phonon-roton gas in the mixture (the nor¬ 
mal 4 He component), the heat capacity of the mixture per mole of 3 He is 
larger for a low 3 He concentration mixture than it is for a high 3 He con¬ 
centration mixture. 14 In the superfluid Stirling cycle, the low concentration 
mixture (with high heat capacity) flows from the low temperature expander 
to the high temperature compressor and the high concentration mixture 
(with low heat capacity) flows from the high temperature compressor to 
the low temperature expander. There is a net cooling effect throughout the 
regenerator due to the imbalance of these heat capacity flows. The tem¬ 
perature distribution in the regenerator is no longer linear. In the extreme 
case where there is no longitudinal heat conduction and the heat transfer 
coefficient between the regenerator and fluid is infinite, the temperature 
profile in the regenerator is uniformly at the expander temperature. Only 
differentially close to the compressor end of the regenerator does the 
regenerator temperature approach the compressor temperature. We con¬ 
clude then that the temperature of the fluid entering and exiting the 
expander end of the ideal regenerator is always at the expander tem¬ 
perature. 

The corresponding comment for the compressor end of the regenerator 
is not true. There is a significant temperature difference between the fluid 
entering and exiting the compressor due to the large imbalance in the heat 
capacities. In the PR model, it is the cold fluid stream exiting the 
regenerator into the hot compressor piston and the hot fluid stream exiting 
the hot compressor piston into the cold regenerator that produces the 
irreversibilities of the cycle. 

The imbalance in the heat capacities can be understood by considering 
the normal part of the 4 He component of the mixture, the phonons and the 
rotons, that are swept along with the 3 He. 15 The phonon-roton gas 
volumetric density is a function only of the temperature and not the 3 He 
concentration. When 3 He flows from the compressor to the expander at 
high concentration, the 3 He sweeps a fairly small volume of phonons and 
rotons through the recuperator and therefore, relatively few phonons and 
rotons need be annihilated to cool the mixture from the compressor to 
expander temperature. On the other hand, when the 3 He flows from the 
expander to the compressor at low concentration, a large volume of 
phonons and rotons are swept through the recuperator and therefore, 
many phonons and rotons are created in warming the mixture from the 
expander to the compressor temperature. Since the energy and entropy 
absorbed (rejected) by the 4 He component of the mixture is proportional 
to the number of created (annihilated) phonons and rotons, we see then 
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that the low concentration flow has a higher heat capacity per 3 He atom 
than the high concentration flow. 

In a real system, the temperature distributions in the regenerator are 
smoother than described in the ideal case due to imperfect transverse con¬ 
duction to the regenerator walls and to longitudinal conduction. 


Cooling Power 

With the temperature of the flows entering and exiting the expander 
piston determined to be T e , the two remaining cyclic integrals can be 
calculated and the cooling power of the phonon-roton model can be 
evaluated. In the ideal gas model, the enthalpy of the 3 He quasiparticle gas, 
h q , is just 5/2 RT. Incorporating the observations of the past few sections 
into Eq. (39) we find 

Qe = W e - 5 -RT e §dN 3 - 7>°<f ^ dN 3 =W e - T e stf^ dN 3 . (43) 


The second equality is true because § dN 3 = 0 by conservation of mass. The 
last integral in Eq. (43) can be expanded as 



( Ve dn 3 T n 3 uf V e ) 



( V e dn s + « 3 dV e ) 


dv, 

„0 


V e dn , 




- — ot 



(44) 


In the final expression of Eq. (44), the first integral is zero because the final 
volume is the same as the initial volume. The parenthetic sum multiplying 
a' is zero by conservation of the 3 He mass. Using these results Eq. (43) 
becomes 


Q .-w.-?gf v - dn > 


vi J n. 


(45) 


The cyclic integral simplifies to 

V ‘ C - 3 = nb V swe sin(/? - a) 


b 2 \yr^b 


+ ndv swe sin(a) ( 


Vi ~ d2 d V* d 


- — 1 


(46) 
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where 


A P f V swc 
2T c (NR) efAe ' 

The total cooling power of the refrigerator is easily calculated by substitut¬ 
ing Eqs. (37) and (46) into Eq. (45), 


Heat Transfer in the Compressor 

The heat transfer experienced by the compressor is most easily 
evaluated by applying the conservation of energy to the entire SSR, 

Qc= ~Qe+W c +W e (47) 

where —Q c is the heat rejected by the SSR compressor to the environment. 
All the quantities on the right hand sides of the Eq. (47) are easily deter¬ 
mined using Eqs. (36), (37), (45) and (46). The COP of the refrigerator is 
determined using the definition contained in Eq. (12). Finally, we define a 
second law efficiency as the ratio of the refrigerator COP to that of a 
reversible refrigerator operating between the same two temperatures: 


> 1 2 = 


COP 

COP„, 


= COP 



(48) 


Entropy and Fountain Pressures of Pure 4 He 

We have now generated expressions for the compressor and expander 
works and heat transfers. These expressions depend on the volumes and 
temperatures of the compressor and expander, the number of moles of 3 He, 
the fountain pressure difference between the expander and compressor, and 
the entropy of pure 4 He at the expander temperature (s® in Eq. (45)). Once 
again, we are interested in developing a closed form expressions for the 
work and heat transfers of the SSR. To this end then, it is necessary to 
have closed form expressions for the entropy and fountain pressures of pure 
4 He. 

We have fit the entropy as a function of temperature to published 
entropy values of pure 4 EIe 16 with an equation of the form: 



e~^ T . 


= 4BT 3 + ^ (y + 5 ) ^ ^ A,/T + C | 


(49) 
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The fit parameters A, B. C, A, and A 2 have values of 23.2 J/mole K 3/2 , 
6.75 x 10 -3 J/mole K 4 , 500 J/mole, 8.65 K and 15.7 K, respectively. This 
expression is a modified form of the entropy expression given by Khalat- 
nikov 17 for the phonon-roton gas present in pure 4 He. The first term is the 
contribution of phonons to the 4 He entropy. The second term is the con¬ 
tribution of the rotons to the entropy of the 4 He, A, is the characteristic 
energy associated with a roton excitation. The last term can be thought of 
as the contribution of the “maxons” to the entropy of the gas and is 
necessary to improve the high temperature, 1.2 K < T < 2 K, fit of the data. 
A graph of the fit is shown in Fig. 3. The relative error between the 
Eq. (49) and the values given by Wilkes is typically less than 4% and never 
exceeds 8 %. 

The data used in determining the coefficients used in Eq. (49) are 
entropies along the saturated vapor curve. We have assumed that the 
entropy is not a function of pressure throughout the range of interest. For 
liquid 4 He at temperatures between 1 and 2 K, a pressure change from one 
to two atmospheres induces a change in the entropy of the fluid of less than 
1%. 18 The pressure terms that were neglected in formulating Eq. (49) are 
smaller than the error associated with the fit. The entropy-independence-of- 
pressure assumption is also consistent with our analysis of the SSR where 



Fig. 3. A plot of the molar entropy of pure 4 He as a function of temperature. The 
dark triangles are experimentally derived values taken from a table of molar 
entropy in Wilkes. The crosses connected by lines are the result of the fit of Eq. (49) 
to these triangle points. 
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we have assumed that the helium mixture behaves as an incompressible 
fluid. 

The fountain pressure of pure 4 He is linked to the entropy relation by 
London’s expression: 


Pf(T)=\ T pls'dT (50) 

J o 

which was used in Eqs. (17) and (18) as the definition of fountain pressure. 
Using Eqs. (49) and (50) and assuming the density of the liquid is inde¬ 
pendent of temperature, the fountain pressure is 


P/ {T) = p° 4 {BT 4 + + Ce ( -^ /T) ). (51) 

With these expressions for the fountain pressure and the entropy, we now 
have a set of analytic expressions for the work and heat transfers of an 
SSR. 


Experimental Apparatus 

The experimental apparatus used in this work is similar to that dis¬ 
cussed in detail in Patel and Brisson 19 so that we will only briefly describe 
the apparatus here. A drawing of the apparatus is shown in Fig. 4. The 
compressor is cooled using a 4 He evaporation refrigerator (1 K pot) which 
is not shown in the figure. 

This SSR is of the Brisson and Swift type having two SSR’s operating 
180° out of phase with each other and using a recuperator. The compressor 
and expander piston platforms of this SSR are solid blocks of OFHC cop¬ 
per on which pistons of edge welded stainless steel bellows are mounted. 
The effective areas and strokes of the compressor pistons are 17.74 cm 2 and 
1.0 cm, respectively, and for the expander pistons 13.61 cm 2 and 0.98 cm. 
Within each piston platform, there are isothermal heat exchangers made 
from nested OFHC copper cylinders press fit into the piston platforms. 
A 76j«m gap exists between the inner wall of an outer cylinder and the 
outer wall of the inner cylinder. In total, there is 65.94 cm 2 of heat transfer 
area in the compressor piston platform and 209.70 cm 2 of heat transfer area 
in the expander piston platform. The recuperator of this SSR is made from 
Kapton and consists of two-hundred 2.38 mm by 127 nm by 20 cm long rec¬ 
tangular passages arranged as a counterflow heat exchanger and separated 
by a 25 ^m Kapton septum. The total volume of the recuperator is 20.18 cm 3 
of which 12.10 cm 3 (6.05 cm 3 per side) is devoted to recuperative heat trans¬ 
fer. The total volume of the SSR is 87.2 cm 3 (43.6 cm 3 per side). 
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Fig, 4. A drawing of the superfluid Stirling refrigerator used in this investigation. 
The pistons are made with bellows to avoid the leaks and the friction associated 
with sliding seals. The recuperator is a plastic composite and the isothermal heat 
exchangers are made from press fit concentric copper cylinders. A pumped 4 He 
refrigerator (not shown) provides cooling for the compressor pistons. 


The SSR was operated with a cycle period of 20 seconds and with a 
3 % 3 He- 4 He mixture. Data was taken for several compressor temperatures. 
The gross control of the compressor temperature was set by throttling the 
4 He flow from the 1 K pot. The fine temperature control was done using 
a feedback controlled electric heater on the compressor platform. Fixed 
heat loads were placed on the expander platform using an electric heater. 
For each heat load, the expander platform temperature was allowed to 
come to steady state and then the expander platform temperature was 
recorded. The compressor and expander temperatures were determined 
using calibrated ruthenium oxide resistors. 20 
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Model versus Experiment 

There is some ambiguity as to what volume values should be chosen 
to model our experimental SSR. The superleaks themselves are 28% void 
space and hence can contain significant amounts of mixture. However, we 
assume that the 3 He component that has diffused into the superleak is 
immobile on the time scale of a single Stirling oscillation. We also assume 
that the 3 He- 4 He mixture concentration in the Vycor superleak is the same 
as that of the bulk mixture in the SSR. Under these assumptions, the effect 
of the superleak on the SSR performance can be neglected. Practically 
speaking, the Vycor superleaks require the loading of some extra mixture 
into the SSR to fill them but they do not alter the concentration or the 
dynamics of the 3 He component of the mixture. 

For one side of the SSR, the compressor clearance and swept volumes 
are 12.04 and 17.74 cm 3 , respectively. Similarly, the expander clearance and 
swept volumes are 9.83 and 13.34 cm 3 , respectively. The clearance volumes 
include the isothermal heat exchanger volumes and any “connecting” 
(headering) volumes between the piston, the isothermal heat exchanger and 
the recuperator heat exchange surfaces. The volume associated with the 
recuperator heat exchange surfaces (6.05 cm 3 ) can be distributed between 
the clearance volumes of the compressor and expander pistons. In the 
calculations that follow, we have arbitrarily added half of this volume to 
the compressor clearance volume and half to the expander clearance volume. 

All the calculated cooling powers and works reported in this article 
are double what would be calculated from the Schmidt and PR model 
expressions. This was done to account for the two back-to-back SSR’s used 
in the experiment. 


Comparison of the PR Model to the Schmidt Model and Experiment 

Figure 5 shows plots of experimental cooling powers, Schmidt model 
cooling power and the PR model cooling power for the superfluid Stirling 
refrigerator for a compressor operating temperature of 1.2 K. In this case, 
the fountain pressure is significantly smaller than the osmotic pressure 
associated with the 3 He component of the mixture so that the influence of 
the phonon-roton gas on the motion of the 3 He gas is minimal. The 
reduced cooling power of the PR model compared to that of the Schmidt 
model at temperatures below 0.8 K is primarily due to the reduced mass 
flow of 3 He between the compressor and expander stages in the PR model. 
(The phonon-roton gas in the compressor piston-cylinder reduces the den¬ 
sity of the 3 He in the compressor and increases the density in the expander 
piston.) 
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Fig. 5. A plot of the SSR cooling power versus expander temperature. The compressor 
is operated at 1.2 K. The three traces correspond to experiment (circles), the Schmidt 
model (squares) and PR model (*). A 3% 3 He in 4 He mixture is assumed. 


The real performance of Stirling refrigerators is always significantly 
worse than the performance predicted by the Schmidt model, since this 
model does not consider any loss mechanisms. The high temperature per¬ 
formance of the SSR in Fig. 5 is a rare example of a real Stirling machine 
outperforming the predictions of the Schmidt model and is a clear indicator 
that the ideal gas model is an inadequate representation of the high tem¬ 
perature behavior of the working fluid. The large increase seen in both the 
experimental and the PR model cooling power curves at high temperature 
comes from the additional entropy flow due to the phonon-roton gas. 
When the expander is at high temperature (f t >0.9K), significant num¬ 
bers of phonon and rotons are generated in the expansion process. These 
phonons and rotons are subsequently swept up the recuperator with the 
3 He quasiparticle gas to the compressor. These phonons and rotons are 
then destroyed in the compression process, releasing their energy and 
entropy to the high temperature heat exchangers. 

It should be noted that the heightened cooling power of the PR model 
over that of the Schmidt model does not constitute a violation of the 
second law because there is a net increase in the work required to drive the 
SSR. We will show that the presence of the phonon-roton gas serves at 
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high temperature to increase the cooling power of the SSR at a cost of an 
increase in entropy production. 

We attribute the large differences at low temperature between the experi¬ 
ment and the PR model to several effects. The first is that the thermal conduc¬ 
tance (the UA product) between the fluid streams in the recuperator is not 
infinite. There are losses in the real machine due to the ineffectiveness of the 
heat exchangers. We would expect these losses to scale with the difference of 
the compressor and expander temperatures. The experimental data agrees 
well with the PR model when this temperature difference is small (0.8 K < 
T e < 1.2 K in Fig. 5). Second, due to Fermi effects, the equation of state of the 
3 He component of the mixture will deviate from the Boltzmann ideal gas law 
at low temperature. Third, the actual expansions and compressions in the 
piston-cylinders are not entirely isothermal as assumed in the PR model. 

The error bars shown on the PR curve in Fig. 5 (and in Figs. 6-8) 
correspond to the variation of the calculated cooling power when the 
recuperator volume is incorporated entirely into the compressor clearance 
volume in one extreme and entirely into the expander clearance volume in 
the other extreme. 

Figure 6 consists of plots of experimental cooling power, Schmidt 
model cooling power and the PR model cooling power for the superfluid 



Fig. 6. A plot of the SSR cooling power versus expander temperature. The compressor 
is operated at 1.4 K, The three traces correspond to experiment (circles), the Schmidt 
model (squares) and PR model (*). A 3% 'He in 4 He mixture is assumed. 
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Stirling refrigerator for a compressor operating temperature of 1.4 K. As 
expected, the Schmidt model cooling power is linear with the expander tem¬ 
perature. The PR model predicts a lower cooling power than the Schmidt 
model due to the flushing of some of the 3 He out of the high temperature 
piston and the reduction of the ’He mass flow through the recuperator. 
Notice that this effect is more pronounced with a compressor temperature 
of 1.4 K (Fig. 6) than with a compressor temperature of 1.2 K (Fig. 5). Con¬ 
trary to the Schmidt model, the PR model predicts that there is a tempera¬ 
ture at which the cooling power of the refrigerator is zero. Extrapolating the 
PR trace in Fig. 6 results in an estimated minimum operating temperature 
for the refrigerator is below 0.2 K as opposed to the experimental data that 
shows a minimum temperature of 0.621 K. We will discuss the minimum 
operating temperature of the SSR in more depth shortly. 

Figure 7 consists of plots of experimental cooling power, Schmidt 
model cooling power, and the PR model cooling power for the superfluid 
Stirling refrigerator for a compressor operating temperature of 1.6 K. The 



Fig. 7. A plot of the SSR cooling power versus expander temperature. The compressor is 
operated at 1.6 K. The three traces correspond to experiment (circles), the Schmidt model 
(squares) and PR model (*). A 3% 3 He in 4 He mixture is assumed. 
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deviation of the PR model from the Schmidt model is more pronounced 
than in the two earlier cases. The PR model minimum operating tem¬ 
perature has increased to 0.5 K versus the experimental minimum 
operating temperature of 1.07 K. 

Finally, Fig. 8 shows plots of experimental cooling power, Schmidt 
model cooling power and the PR model cooling power for the superfluid 
Stirling refrigerator for a compressor operating temperature of 1.8 K. The 
PR model and experimental cooling power curves no longer resemble the 
Schmidt model trace. Both of the PR model and the experimental traces 
show a clear minimum operating temperature of 1.06 K and 1.46 K, respec¬ 
tively. 

The expander and compressor works ( W e , W respectively), cooling 
power (QJ, waste heat (Q,.) and the second law efficiency (f} 2 ) for the PR 
model of the experimental SSR above are shown in Fig. 9. The compressor 
temperature is assumed fixed at 1.8 K and the •'He density corresponds to 
a 3% mixture. 

One of the very clear differences between this model and the Schmidt 
model is the second law efficiency. In the Schmidt model, the second law 
efficiency has a constant value of one; however, in the PR model, the 



Fig. 8. A plot of the SSR cooling power versus expander temperature. The compressor 
is operated at 1.8 K. The three traces correspond to experiment (circles), the Schmidt 
model (squares) and PR model (*). A 3% -‘He in 4 He mixture is assumed. 
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Expander temperature (K) 

Fig. 9. Plots of the expander and compressor works ( W e and W,., respectively), 
expander and compressor heats (Q t , and Q,., respectively) and the second law 
efficiency, as a function of expander temperature for the PR model. The 

compressor temperature is assumed to be 1.8 K and the loaded 3 He concentration 
is 3%. 


second law efficiency drops to zero due to the heat capacity flow mismatch 
in the recuperator. Figure 10 shows the second law efficiency of the 
modeled SSR as a function of expander temperature for several compressor 
temperatures. On this graph, a Schmidt model second law efficiency would 
appear as a horizontal line with a value of r] 2 - 1. As expected, in the low 
temperature compressor limit ( T c = 0.8 K), the PR model efficiency 
approaches that of the Schmidt model. The plots of efficiency become 
suspect at very low temperatures 7^ 0.5 K due to the change of the 3 He 
from Boltzmann gas behavior to Fermi gas behavior. 


Ultimate Temperatures 

The minimum expander temperature achievable by an SSR is a func¬ 
tion of the compressor temperature, the average 3 He concentration and the 
ratio of the compressor and expander volumes. Practically speaking, the 
ultimate temperature is always influenced by the recuperator effectiveness 
and expander efficiency; but here, we are interested in determining an 
ultimate temperature based on 100% effective components. 
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Expander temperature (K) 

Fig. 10. Parametric curves for the second law efficiencies for the PR model as a function of 
the expander temperature. The volumes assumed correspond to the experimental apparatus. 


The minimum achievable temperature is linked to the flushing of all 
the 3 He from the compressor piston into the expander piston due to the 
fountain effect. When this happens, there is a column of high-thermal-con¬ 
ductivity pure 4 He between the high temperature compressor and the low 
temperature expander which thermally “shorts” the refrigerator. The 
expander temperature at which this occurs can be determined by setting p qc 
in Eq. (21) equal to zero. The result is a transcendental expression for 
T e | min , the minimum expander operating temperature: 


(p f (T c )~p f (T e LJ) 
N,R 


(52) 


where V e | max is the maximum volume of the expander ( = V ch , + V sn , e ). In 
general, the SSR is loaded with a mixture of known concentration when 
both the compressor and expander pistons are in their center positions and 
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at uniform temperature. This means that total number of moles of 3 He can 
be written as: 


N* = n-. 


v clc +^+v cle +-f 


Vclc + + Vde + ~^) X load 


o2(l + 


where n 3 load and x, oad are the molar 3 He density and the molar concentra¬ 
tion of the mixture when it is loaded into the SSR. Substituting Eq. (53) 
into Eq. (52) we find: 


_ i f T t yni _ _ (1 + x'x /oad ) l 

. Vetc + ( f^swc/2) + V cle + ( V swe /2) x /oa j 

: . v° 4 (Pf(T c )-p f (T e 1^)) 

R 

g A(PA T c)~Pf(T e Imin)) 

^ r> 


where the dimensionless quantity £ is defined as the quantity in the square 
brackets. is a product of two terms. The first is the ratio of maximum 
expander volume to the total (loading) volume of the SSR and is deter¬ 
mined entirely by the geometry of the SSR. The second term is (to first 
order) the reciprocal of the molar concentration of 3 He loaded into the 
SSR. In the SSR’s built to date, the values of £ has typically varied from 
1 to 100. Solutions for T e |„„„ using Eqs. (51) and (54) are plotted in Fig. 11 
in terms of £ and for T c values of 2.0, 1.8, 1.6, 1.4, 1.2, and 1.0 K. 

Figure 11 is a quick reference to estimate the ultimate temperature of 
a superfluid Stirling machine given the overall geometry, the 3 He concen¬ 
tration loaded into the refrigerator, and the compressor operating tem¬ 
perature. For example, the SSR used in the experimental part of this work 
was loaded with a 3% mixture and has a maximum expander volume to 
loading volume ratio of 0.603 which results in a £ value of 20.2. Using 
Fig. 11, we find that the ultimate operating temperature is limited to 1.72 K 
for 7 C = 2K, to 1.31 K for T c = 1.8 K, to 0.68 K for T c = 1.6 K, to 0.27 K 
for r c =1.4K, to 0.09 K for T c = 1.2 K, and to 0.027 K for T C =1.0K. 
These ultimate temperatures are indicated in Figs. 5 through 8 by a dark 
arrow. A more conservative estimate can be made by including all the 
recuperator volume as part of the maximum expander volume. Then, the 
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Fig. 11. A plot of the theoretically achievable expander temperature due to phonon- 
roton effects as a function of the dimensionless parameter ( using Eq. (54) (solid lines). 
Each curve is labeled with the compressor operating temperature. The A-line (dotted), 
the dilute phase line (dash-dot), and the ’He Fermi temperature (dash-dot-dot) are 
also plotted on the same axes. The labeled vertical dashed lines indicate the mixture 
concentration in the expander for a given ultimate expander temperature and q. 


maximum expander volume to loading volume ratio becomes 0.672 and £ 
has a value of 22.6. Using Fig. 11, we find that the ultimate operating tem¬ 
perature is limited to 1.75 K for T, = 2 K, 1.36 K for T, = 1.8 K, to 0.76 K 
for T,.= 1.6 K, to 0.31 K for T,.= 1.4 K, to 0.10 K for T,.= 1.2 K, and to 
0.03 K for T,.= 1.0 K. These ultimate temperatures are indicated in Figs. 5 
through 8 by a light arrow. There is an apparent contradiction in Figs. 5 
through 8 between the ultimate temperature calculated using Eq. (54) 
(Fig. 11) and the calculated zero cooling power using Eq. (45). In deriving 
Eq. (45) no limits were imposed on the sign of the pressure and number 
density of the 3 He. Under certain operating conditions, it is possible (in the 
PR model, not the experiment) to have a negative mass density of 3 He. 
This is clearly unphysical and the solutions of Eq. (54) determine the point 
when the PR model breaks down. 
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It becomes clear from Fig. 11 that high compressor temperatures 
necessitate operating an SSR with a small value of £ to achieve low tem¬ 
perature. To achieve small values of £, the volume of the expander must be 
small fraction of the total volume of the SSR and the loading concentration 
of the mixture needs to be high. 

In our SSR, the heat flush effects are largely irrelevant for determining 
the ultimate temperature for compressor temperatures between 1.2-1.6K, 
as the low temperature performance of the SSR is dominated by the 
recuperator ineffectiveness. The heat flush effect is clearly seen in Fig. 8 
where the compressor temperature is 1.8 K. The experimental curve is quite 
smooth until it achieves a temperature of 1.46 K where there is a definite 
“knee” in the cooling power curve as it drops abruptly to zero cooling 
power. Strictly speaking then, the PR model curves shown in Figs. 5-8 
each should have knees where the cooling power abruptly drops to zero at 
the “arrowed” temperature. We believe that the difference between the knee 
temperature and the arrowed temperature in Fig. 8 can be attributed to the 
inaccuracies in the fountain pressure fit, in the values of the SSR volumes, 
and in the concentration of the mixture loaded into the SSR. 

A consequence of a large recuperator volume is the broadening of the 
“heat flush transition” in the Q c vs. T e plot. As the expander temperature 
decreases, the region of pure superfluid 4 He begins to oscillate into and out 
of the recuperator, decreasing its effectiveness. If the recuperator has a large 
volume, there will be a large expander temperature range over which the 
thermal isolation afforded by the recuperator drops to zero. 

The operation of the SSR is limited to the superfluid range of the 
phase diagram to allow superflow through the superleaks. The two-phase 
region is forbidden. As the SSR cools, the concentration of the 3 He 
increases in the expander. If the temperature is low enough and the concen¬ 
tration is high enough the mixture will spontaneously phase separate in the 
expander. This leads to a “locking” of the dilute 3 He concentration in the 
expander to the concentration associated with the phase transition at the 
expander temperature, precluding further cooling with the Stirling cycle. 

Estimates on these limits can be accounted for in Fig. 11. ^ can be 
interpreted as a non-dimensional 3 He molar density at the presumed 
ultimate operating point. Using Eqs. (53) and (41), £ can be written as 

£ _ Ftejmax (55) 

V° 

where v 3e | max is the molar volume of the 3 He when all the 3 He is placed 
into the maximum volume of the expander piston. £ is then a comparison 
of the molar volume of the 3 He, when all of the 3 He in the SSR is in the 
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expander, to the molar volume of pure 4 He. Using Eq. (55), the superfluid 
transition for the mixture (the A-line) and the dilute-phase two-phase 
boundary 21 are also plotted in Fig. 11. SSR’s designed with £’s to the left 
of these A and two phase curves will have ultimate temperatures above that 
predicted by the solid curves in Fig. 11 because these (£, T) points are inac¬ 
cessible to an operating SSR. The ultimate operating temperature of these 
machines can be determined by finding the maximum 3 He quasiparticle 
density in the expander as a function of expander temperature, using Eq. 
(26) and the ideal gas law. Where this curve crosses the dilute two phase 
line determines the estimate for the ultimate temperature of that design. 

We have also plotted the Fermi temperature of the 3 He as a function 
of in Fig. 11 as a double-dot single-dash curve. This curve is an indication 
of where the 3 He Fermi equation of state begins to deviate from that of an 
ideal Boltzmann gas and hence, marks where the PR model begins to 
break down. For our experimental SSR, the ultimate temperature estimates 
made earlier for compressor temperatures less than 1.4 K are rather suspect 
in that all the £,, T e 1^ points fall well below and to the left of the Fermi 
temperature line. In this region, the Fermi gas model is more appropriate 
than the Boltzmann gas model for the 3 He quasiparticle gas. In generating 
the Fermi temperature curve we have assumed that the effective mass of the 
3 EIe atoms in solution are 2.5 times the mass of a bare 3 He atom. 22 Finally, 
we have included four vertical dashed lines, labeled 20%, 10%, 5% and 
1% to indicate the 3 He concentration in the expander when an SSR (of a 
given <J) is operated at its 7 , <f | min . 


Closing Comments 

The performance limitations on the SSR described here are ther¬ 
modynamic, imposed by the non-ideal thermodynamic properties of the 
working fluid. These limitations are not seen with the Schmidt model, 
which predicts an ultimate expander temperature of zero Kelvin. All 
assumptions concerning the transport properties have been ideal. For 
example, we have implicitly assumed that the viscosity is zero, that the 
thermal conductivity in the recuperator is infinite in the direction normal 
to the flow and zero in the direction parallel to the flow, and that the ther¬ 
mal contact in the piston cylinders to the thermal reservoirs is perfect. The 
overall performance of the real SSR will be worse than what is predicted 
by the PR model depending on the transport losses. 

The PR model is simple and the equations of state that we have used 
above are not appropriate over the entire possible operating range of the 
SSR. At high concentrations, for example the 3 He does not behave entirely 
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as an ideal gas. In addition, the 4 He properties change as the 3 He concen¬ 
tration is increased. SSR’s have been operated with 3 He concentrations as 
high as 36%, 23,24 but, it is very questionable whether the ideal gas model 
is appropriate for these concentrations. Nevertheless, the PR model does 
very well in capturing many of the features and magnitudes of SSR perfor¬ 
mance for temperatures between 0.5 K and 2 K. This model is a useful tool 
to design an SSR and provides a useful handle on performance limitations. 

Once again, we have outlined a new simple closed-form model for the 
high temperature operation of the superfluid Stirling refrigerator. This 
model is intended as a design tool to estimate the cooling power and 
ultimate temperature of any given design. The model shows that for tem¬ 
peratures above one Kelvin, the 4 He component significantly alters the per¬ 
formance of the SSR from that of standard Stirling machines. The National 
Science Foundation and the McCarthy chair supported this work. 
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